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We investigate coexistence of antiferromagnetic and superconducting correlations in bilayered ma- 
terials using a two-dimensional t-J model with couplings across the layers using variational Monte 
Carlo calculations. It is found that the underdoped regime supports a coexisting phase, beyond 
which the (d-wave) superconducting state becomes stable. Further, the effects of interplanar cou- 
pling parameters on the magnetic and superconducting correlations as a function of hole doping 
are studied in details. The magnetic correlations are found to diminish with increasing interplanar 
hopping away from half filling, while the exchange across the layers strengthens interplanar antifer- 
romagnetic correlations both at and away from half filling. The superconducting correlations show 
more interesting features where larger interplanar hopping considerably reduces planar correlations 
at optimal doping, while an opposite behaviour, i.e. stabilisation of the superconducting state is 
realised in the overdoped regime, with the interplanar exchange all the while playing a dormant 
role. 
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I. INTRODUCTION 

The mechanism leading to electron pairing in cop- 
per oxide superconductors has stimulated a great deal 
of speculation. The manifestation of planar antiferro- 
magnetism in the Cu02 layers have provided motivation 
to study the role and importance of antiferromagnetic 
(AF) interactions between Cu^^ spins and their intimate 
relevance to superconductivity^. Doping with holes in 
these insulating cuprates results in destruction of long 
range order, however short range antiferromagnetic cor- 
relations between the copper moments survive. Further 
doping leads to the emergence of a superconducting (SC) 
state. 

Thus the interplay of AF and SC phases has gener- 
ated much attention. A phenomenological 5*0(5) the- 
ory has attempted to unify antiferromagnetism and su- 
perconductivity owing to their proximity in the phase 
diagrams.!-. The basic assumption of the theory is that 
these two phases share a common microscopic origin and 
hence both demand treatment at an equal footing. A 
number of experimentally observed features, such as the 
vortex state, a resonant peak in SC state from neutron 
scattering data for optimally and underdoped samples 
etc. have provided ample credence to the theory. How- 
ever a microscopic theory in this regard is still lacking. 

A natural extrapolation of finding the connection be- 
tween AF and SC is the issue of their coexistence which 
seems more crucial and calls for attention. A large vol- 
ume of work exists that focuses on the various details 
of the coexistence phenomenon. Some of the theoretical 
attempts include mean field studies of t-J and Hubbard- 
like models which confirm the coexistence of magnetic 
and superconducting orderi^. The coexistence is sug- 
gestive of the presence of short range AF correlations in 
SC state that are probed by inelastic neutron scattering 
experiments via an enhanced scattering intensity near the 



AF wavevector (tt, 7r)^iii^i^. 

The coexistence issue has been revived recently in the 
context of bilayer (and multilayer) cuprates^. While 
acknowledging the planar correlations perhaps dominate 
the physical properties of these superconductors, the role 
of interlayer couplings and their relevance to the coexis- 
tence phenomena for layered materials have been heavily 
emphasized. Intimately connected to this is the question: 
whether the superconducting correlations originate from 
AF spin fluctuationsi2iiiiii^ or via a electron-phonon me- 
diated pairing enhanced by interlayer tunneling^. How- 
ever both these mechanisms cannot be operative to- 
gether. 

The next fundamental question is the symmetry of the 
superconducting gap function. Even with sufficient ex- 
perimental evidence for a dx^^-y^ pairing for a planar ma- 
terials (for a comprehensive review on the subject see 
Ref. [3), for bilayers the pairing symmetry is still un- 
clear. Several variational calculation o^^i^^'^^ performed 
over the years have provided nourishment to a dx2_y2- 
wave pairing scenario and a reasonably broad window 
of carrier concentration has been identified in planar sys- 
tems where magnetic and SC order coexist. Similar stud- 
ies in the context of bilayers are lacking and hence provide 
motivation for us to investigate a bilayer t-J model via 
variational Monte Carlo (VMC) technique. 

Our goal in this work is to examine the coexistence of 
antiferromagnetism and superconductivity in t-J bilayers 
using VMC and to study the magnetic and superconduct- 
ing correlation in the variational ground state. We fur- 
ther intend to investigate the dependence of these proper- 
ties on interlayer coupling strengths. In a recent paper—, 
to determine the most suitable pairing symmetry of the 
SC state in t-J bilayers we investigated the stability of 
various pairing symmetries, e.g. (a) A(cosfca; — cosfcj,) 
(d-wave), (b) A|| (cosfcj; — cosfcy) -|- Aj_ cosfcj and (c) 
A| I (cosfcj, — cosfcj,) 4- A^(l — cosfcz). Another function 
discussed in connection with bilayer materials but not 
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included in our previous study is the state (A^ — 
±A(cosfc2: + cosfcy) with '+' sign for — and '— ' 
for fez = tt) which has s-wave symmetry and opposite 
signs in the bonding (fe^ = 0) and antibonding bands 
{kz = 7r)^ii^. This state possesses the merit of explain- 
ing the resonance peak at 41 meV obtained in neutron 
scattering experiments^i^Si^iiS^. In this paper we con- 
sider also this pairing symmetry in the search for most 
stable ground state in bilayers. 

Here we state the main results obtained by us. The 
long range AF order coexists with superconductivity in 
the underdoped regime. In the coexisting phase not 
only the AF but SC correlations are also significantly 
stronger as compared to that in the pure SC (d-wave) 
state. Larger interlayer hopping frequency reduces pla- 
nar SC correlations in the optimally doped phase whereas 
it enhances it in the overdoped regime, while the effect of 
interlayer exchange on SC correlations is minimal for the 
range of the parameter values considered in our paper. 

Our paper is organised as follows: section II introduces 
the t-J model for bilayers and discusses the most suitable 
variational wavefunction to be used for our calculations. 
A brief note on the numerics used and an elaborate dis- 
cussion on the results appear in section III. The effects 
of interplanar coupling parameters have been emphasized 
in magnetic, superconducting and the coexisting phase. 
Section IV concludes with a brief summary of the results 
obtained in this paper. 

II. THE HAMILTONIAN AND THE 
VARIATIONAL WAVEFUNCTION 



The t-J Hamiltonian for a bilayer can be written as 




where t and J are the planar hopping and exchange inte- 
gral respectively, while t± and J± are the corresponding 
interplanar parameters. Qo- (c|^) annihilates (creates) 
an electron of spin a at site i, rii = cla'^icr a-nd S; 
is the spin operator at site i given by Sf — ipj{^cra)ipi. 
ipl = ^c|| and cTq, (with a — x,y, z) are the Pauli 

spin matrices. The summation indices and {{i,k)) 
indicate nearest neighbour pairs in the same plane and 
different planes respectively. The Hamiltonian obeys an 
essential requirement, i.e. it acts on a subspace of no 
doubly occupied sites. 

To incorporate the coexistence of AF and SC phases 
we consider the following variational wavefunction as the 
ground state of the Hamiltonian. 

l^yar (Asc, Aaf)) = VgVn J] ("k + ^'k^T^'-ki) |0> (2) 

k 



where the operator Vn projects out the states with a 
fixed electron number, N and Vg — ni(l~'^iT'^ii)' 
Gutzwiller projector which imposes the condition of no 
double occupancy. The product in Eq.[2]is for over all the 
'k' points in the first Brillouin zone and the amplitudes 
Uk and Vk are defined by, 

^ = ^(k) = (3) 

(Ti^k - A^) + V(T£;k - /^)' + Al 

where Ak = Asc/(k) represents the SC gap, /(k) be- 
ing an appropriate symmetry function of k and E]^ = 

y^e| + A^y. ek = — 2i (cos kx + cos ky ) — 2t± cos kz , is the 
free electron dispersion and /i is the chemical potential. 
The —(-!-) signs in the denominator of Eq. [3] corresponds 
to Ek < (ek > 0). The quasiparticle operators, dj^^ 
diagonalizes the AF Hartree-Fock Hamiltoniartii with a 
gap Aaf and are related to the electron operators by the 
following transformation. 




Here Q = (tt, tt, tt), is the perfect nesting vector and 

r/(a) = ±lforfT=r,i- 

The wavefunction in Eq. [2] consists of two variational 
parameters, viz. Age and Aaf. Ideally the chemical po- 
tential, /i should also be treated as a variational parame- 
ter, however here we fix it at its noninteracting value, fiQ. 
This is because the energy correction obtained by vary- 
ing /i has been found to be negligibly small (for small 
J)^^ for a square lattice and we expect it to be the same 
also for bilayers. The wavefunction describes different 
phases depending upon the (relative) values of the vari- 
ational parameters. For example, Aaf = describes the 
usual BCS superconducting state^i, whereas in the limit 
Age 0, the wavefunction reduces to a state with anti- 
ferromagnetic long range order—. For nonzero Age and 
Aaf , the wavefunction describes a phase with coexisting 
AF and SC state, while the normal state is recovered as 
both the parameters vanish. 

III. THE RESULTS 

We shall skip details of the variational Monte Carlo 
method used as it appears elsewher o^'^i^^'^^ and only 
provide the essential features of our computation. We 
consider periodic boundary condition along the planar 
x-direction and antiperiodic boundary condition in the 
planar y-direction to avoid singularity in (/)(k) for the k- 
points with k^ = ky and ei!_ — fJ, < 0^^. In one Monte 
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Carlo Sweep (MCS) through the lattice, Ng (equal to 
the number of lattice sites) random moves are attempted 
which consists of moving an electron to an empty site 
and exchanging two antiparallel spins. After each suc- 
cessful move, Monte Carlo updates of the configurations 
are made by using the inverse update method'^^. Var- 
ious expectation values that are obtained in the paper 
are computed by sampling configurations chosen from 
about 10^ - 10^ MCSs after taking 5,000 warm up sweeps. 
Simulations are performed on a bilayered lattice of size 
10 X 10 X 2. 

In the following, we discuss the stability of different 
phases with different SC pairing symmetries, followed 
by the results for magnetic and superconducting correla- 
tions. We show that the AF and SC phases coexist in the 
underdoped regime by comparing the optimal energy of 
the variational wavefunction having pure SC correlations 
with the one having coexisting AF and SC order. To 
characterize the pure SC and the coexisting phase and 
to enumerate the differences in their properties, we com- 
pute the correlation functions for both these phases and 
make a detailed comparison between the two. 

A. Stability of different phases 

We first consider the pure SC wavefunction (obtained 
by putting = in Eq. [J) . As for the pairing sym- 
metry of the SC state, we considered energies of four 
different variational wavefunctions listed in the previous 
section. An earlier work^'^ investigates in details the first 
three wavefunctions in the list and found that the pure 
d-wave state yields lowest energy at all values of 5 away 
from half filling. This led to the conclusion that of the 
three pairing symmetries discussed there, a planar d- 
wave state is most appropriate in the context of a t-J 
bilayer. We have included the s*-state for comparison 
with the existing ones in the light of the emphasis given 
to it where it is claimed to be more stable than that of 
the d-wave state for the bilayered systems and is capable 
of explaining the origin of the neutron scattering peak 
observed experimentally in YBa2Cu307 at temperatures 
below T^ii2ii2S. However our VMC calculations indicate 
that the energy of the pairing state to be actually 
higher than the normal state at all values of hole concen- 
trations. Thus we discard this pairing symmetry from 
the list of possible candidates and conclude that a planar 
d-wave state is most appropriate to study superconduc- 
tivity in bilayers. 

B. Ground state Energy 

Before we proceed to calculate physical quantities, such 
as the ground state energy etc., it is somewhat interesting 
to look at the variation of the optimal superconducting 
variational parameter, Asc, as a function of hole concen- 
tration for a few representative values of the interplanar 



coupling parameters. The choices of these parameters 
are chosen from experimental dat a'^^i'^^ . From Fig. 1 it 
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FIG. 1: Optimal value of the variational parameter, Asc 
shown as a function of hole concentration, 5 corresponding 
to different choices for the interplanar parameters. All the 
figures in the paper correspond to calculations done on a lat- 
tice size of 10 X 10 X 2. The parameter values are taken in 
units of (planar) hopping t. 

may be noticed that the critical hole concentration, Sc 
at which Asc vanishes, increases (from ~ 0.3 to ~ 0.34) 
with larger interplanar hopping, t±, while the interpla- 
nar exchange, J± has no significant effect on Sc- Below 
Sc, Asc is slightly reduced by both higher values of t± 
and J±. It is worth mentioning here that in a two di- 
mensional square lattice (with same values for the pla- 
nar parameters), a d-wave state is stabilized upto ~ 28% 
hole concentration^^, which is lower than the correspond- 
ing value i.e. 34% obtained here for the bilayer (Fig. 1). 
Thus stability of the superconducting state extends upto 
higher values of hole concentration in bilayers than in 
planar materials. 

Next we introduce the second variational parameter, 
i.e. Aaf into the problem and carry out minimization 
of energy in two-variational parameters space, Asc and 
Aaf. The calculation shows that the energy is signifi- 
cantly lowered in the underdoped regime {S < 0.14) when 
compared to that obtained for the pure SC state. The 
optimal energy, Emin/t (per site) as function of S for the 
two cases is shown in Fig. [2] for one one particular choice 
for the interplanar parameters. The energy difference 
between the two phases is maximum at half-filling and 
decreases gradually with increasing hole concentration, 
finally vanishing at S ^ 0.14. Thus superconductivity 
coexists with antiferromagnetism in the underdoped re- 
gion for a bilayer, a feature also observed for the two 
dimensional t-J model^^. This is one of the key results 
of our paper. Similar energy difference of the two phases 
are found for other choices of t±_ and .J± (as appear in 
Fig. 1). 

C. Magnetic order 

We first examine the magnetic correlations in the pure 
d-wave state. The relevant quantities to compute are pla- 
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FIG. 2: Optimal energy, Emin/t (per site) for the pure d-wave 
SC state and the coexisting AF and SC state as a function of 
hole concentration, 5 for t±^ = 0.20 and J±_ = 0.10. 



nar spin-spin correlations, {S^Sj) and the corresponding 
Fourier transform, S'(q), called as the structure factor 
which is defined as, 



(6) 



The real-space correlations (Fig. 3(a) I shows signature 



of AF order in the planes. However the correlations are 



rable magnitude to that at (7r,7r,7r), which corroborates 
the presence of weak interplanar correlations mentioned 
earlier—. Away from half- filling, both planar and inter- 
planar spin correlations are found to vanish rapidly with 
hole doping. Other values for the interplanar parame- 
ters among the ones considered here are found to have 
little impact on spin correlations except for a small reduc- 
tion of planar correlations by larger interplanar hopping 
away from half filling. The value of staggered structure 
factor, S{'K, TT, tt) at half-filling for the d-wave state is ob- 
tained as ~ 0.03. This value may be compared with the 
exact diagonalization results for a bilayer t-J model^ 
which yields an estimate for S{'k, tt, tt) to be ^ 0.45 for 
Jji_/J~ 0.28 (same corresponding value is used in Fig. ^ 
at half filling^. Thus the staggered magnetization for 
the pure d-wave state obtained in our calculations is far 
lower than that obtained via exact diagonalization stud- 
ies. This large discrepancy can be attributed mainly to 
the absence of interlayer magnetic correlations for the 
pure d-wave state. 

Next we discuss the magnetic correlations in the coex- 
isting AF and d-wave SC state which is lower in energy 
than that of the pure d-wave SC state in the underdoped 
region. In Fig. [4l we show the planar and interplanar 
spin correlations both at and away from half-filling. The 
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(a) Planar spin correlations 
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FIG. 3: (a) Spin-spin correlation as a function of distance 
along a planar edge (x-direction) and along a planar diagonal 
(d-direction) . (b) Spin structure factor, S(q) as a function of 
q chosen along the symmetry lines shown in figure. The plots 
shown are at half-filling and for t±_ = 0.20 and ,J± = 0.10. 

found to decay as a function of planar distances suggest- 
ing an absence of AF long range order (AFLRO) even 
at half filling. It may be noted that the d-wave state 
does not show long range magnetic order also in two di- 
mensional systems. As for the interplanar correlations, 
the d-wave shows AF ordering between the spins in the 
two layers, however the correlations are found to be very 
weak. For instance, the strength of the nearest neighbour 
spin correlation for two corresponding sites in different 
planes is approxima tely 1 0% of that for two sites in the 
same plane. In Fig. |3(b)[ we plot 5(q) as a function of 
{qx, qy), the points being chosen along a symmetry path 
for qz equal to both and tt. The peak in S'(q) at (tt, tt, tt) 
indicates the existence of antiferromagnetic correlations 
in the lattice. However, the peak at (tt, tt, 0) is of compa- 
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FIG. 4: Spin spin correlation as a function of distance along in- 
direction for the coexisting AF and d-wave SC (AF-SC) state. 
The spins are on lattice sites lying (a) along an edge in a plane, 
and (b) on two different planes along the same direction. The 
values of the hole concentration, 5 are shown in the figures. 
Here the interplanar parameters are again chosen as t±_ — 0.20 
and J± = 0.10. 

correlations are clearly much stronger in this case than 
for the pure d-wave state. The interplanar correlations, 
which was very weak in the d-wave state, is almost of the 
same magnitude as the planar correlations. In addition, 
the magnitude of correlations does not seem to decay 
with distance at and even slightly away from half-filling. 
However, at larger values of 6, the magnetic correlations 
decay rapidly with distance as seen for S = 0.14 in Fig.|4l 
The energy calculations also provide a support for this 
result where it is found that the kinetic energy dominates 
over the exchange energy and consequently the antifer- 
romagnetic phase disappears corresponding to i5 '--^ 0.14. 
The presence of magnetic order in the coexisting phase is 
further emphasized by plotting S'(q) versus q at various 
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hole concentrations in Fig.[5l The sharp peaks in S{q) at 
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FIG. 5: S'(q) as a function of q for the AF-SC wavefunction 
at various hole concentrations shown in figure. The {qx,qy) 
points are chosen in the same way as in Fig. 3(b) t± = 0.20 
and J± = 0.10. 



q — (tt, tt, tt) for small values of hole doping (Fig. 5(a) ) 
indicates the existence strong AF long range correlations 
in the system. It should be noted that the value of 
S'(7r, TT, tt) at half-filling is ~ 0.35 vi^hich is considerably 
larger than that obtained for the pure SC state and is 
comparable to the exact diagonalization value (viz. 0.45) 
mentioned earlier—. Thus the support for the coexisting 
phase becomes more robust. Also S'(7r,7r,7r) decreases 
with increasing hole concentration as the magnetic cor- 
relations are weakened by the mobile holes. For qz — 



(Fig. 5(b)), it is observed that S'(q) increases with in- 
creasing hole concentration, the increase being maximum 
at (tt, tt, 0). This signals rapid suppression of interplanar 
AF long range correlations away from half-filling. 

Next we incorporate the effect of interplanar couplings 
on the magnetic correlations. Fig.[6]shows the variations 
of ^(Tr, TT, tt) and ^(Tr, tt, 0) with t± a nd J± both at and 
away from half-filling. In Figs. 6(a) and (b), we let t± 
vary while keeping J± constant. It is seen that t±_ has 
no effect on magnetic correlations at half-filling, which 
is expected as hopping of electron is forbidden at half- 
filling due to the no double occupancy constraint. Away 
from half-filling, S{Tr, tt, tt) shows a decrease with increas- 
ing t± thus indicating AFLRO diminishes by larger in- 
terplanar hopping. Further, the decrease of 5'(7r, 7r,0) 
with increasing tj_ implies that the planar correlations 
are mainly affected. On the other hand, variation of in- 
terplanar exchange is found to have the reverse effect on 
the correlations. Fig 6(c) shows that ^(Tr, tt, tt) increases 
with larger J± both at and away from half-filling. How- 



ever S'(7r, 7r,0) (Fig 6(d) I decreases with increasing J^. 
This is a reflection of the fact that interplanar AF cor- 
relations are strengthened by larger J±, while the pla- 
nar correlations are almost unaffected. This agrees with 
exact diagonalization results obtained for the bilayer t-J 
model^. Thus we conclude that larger t± reduces planar 
magnetic correlations away from half-filling, while larger 
J± enhances interplanar correlations both at and away 
from half-filling. 
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FIG. 6: S{n,n,n) and 5'(7r,7r,0) is plotted as a function of 
hole concentration for various values of interplanar parame- 
ters shown in figure, showing the effect of interplanar cou- 
plings on the AF long range order for the AF-SC state, (a), 
(b) involves variation of t± for a fixed J± and (c), (d) contain 
variation of J± for a fixed t±. 



D. Superconducting correlations 

Our next job constitutes of investigating SC correla- 
tions for both pure SC and coexisting AF and SC state. 
The SC correlation function is defined as, 



(7) 



where B's are pair operators and are represented by 
Br'p — ^(cr'|Cr'+/3| — Cr'xCr'+/3|) which annihilates a 
singlet pair on bond (r',r' -I- /3) and Bl^ creates one 
on (r,r + a), ct and /3 are unit vectors connecting to 
nearest neighbours in the x, y (planar) and z (across the 
plane) directions. We have computed i^a,/3 as a func- 
tion of distance |r — r'| for different hole concentrations 
corresponding to several choices for the interplanar pa- 
rameters. A very useful quantity in this connection is 
the SC order order parameter $ which is obtained as^l 



±<I> for large |r — r'|, with the sign being 



+ (— ) for a to be || (_L) to /3 (both a and /3 lie on a single 
layer). 

To examine interplanar SC correlations, we have cal- 
culated i^a,/3(r — r') taking both a, (3 to be along z- 
direction. The values obtained at all hole concentrations 
are very small (smaller than the error bars), i.e. negligi- 
ble in comparison to the planar correlation values. The 
result appears as no surprise as the pure d-wave state 
which is found to most suitably describe pairing symme- 
try for a bilayer, contains no significant interplanar SC 
correlations. 

Next we show that the planar SC correlations are 
stronger in the AF-SC state in the region of hole dop- 
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ing {i.e. underdoped) in which the two phases coex- 
ist in comparison to that in the pure SC state. Fig. [7] 
shows the SC order parameter, $ corresponding to pla- 
nar Fa,f3{r — r') (i.e. {a, 13) along {x,y) as function of 
S for both the pure SC and the AF-SC state. The fig- 




FIG. 7: SC order parameter, $ shown as a function of hole 
concentration, 5 for the AF-SC and pure SC state. Here the 
interplanar parameters are chosen as t± = 0.20 and J± = 
0.10. 

ure clearly shows that SC correlations are stronger in the 
coexisting phase than in the pure SC state. It is very 
interesting to note that the coexisting phase which is en- 
ergetically favourable in the underdoped region, gives rise 
to not only stronger AF correlations but also enhanced 
SC correlations as compared to that in pure SC state. 
This is another important result of our paper. 

Next to estimate the effects of interplanar couplings 
on the SC correlations in a plane, we calculate the SC 
order parameter, $ as a function of hole concentrations 
for different choices of t± and J± . The results are shown 
in Fig. [8] with t± varying from 0.05 to 0.2 and J± = 0.1. 
From the figure one observes that the SC correlations 
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At optimal values of doping, i.e. in the peak region of 
order parameter <&, larger interplanar hopping reduces 
planar SC correlations strongly. However, just the oppo- 
site behaviour is observed in the overdoped region. As 
shown in Fig. [51 for t± = 0.05, the optimally doped re- 
gion is marked by largest $, indicating greater stability 
of the SC state for smaller t±. Whereas the critical hole 
concentration, 6c at which the phase transition from SC 
to normal state takes place, is ~ 0.30 for t± — 0.05, while 
for t± — 0.10 and 0.20, ScS are obtained as 0.33 and 0.34 
respectively. Thus increasing the interlayer hopping in 
a bilayer results in greater stability of the SC state in 
the overdoped region. With regard to the impact due to 
interplanar exchange, J±, our calculations show that J± 
has no perceptible effect (figures not shown here) on the 
SC correlations at least for the values of J± considered 
here. 

The asymmetric behaviour shown by superconducting 
correlations in the optimally doped and the overdoped 
regime as t± is increased from 0.05 to 0.20 constitutes 
another key finding of this paper. While the reduction of 
planar SC correlations with increasing t± as observed in 
optimally doped region is expected due to the fact that 
larger hopping have a disrupting on the pairing, the rea- 
son for the opposite behaviour in the overdoped regime 
where SC correlations are enhanced by larger t±, is not 
immediately obvious. To explain this we look at mo- 
mentum distribution for two different values of t±, viz. 
0.05 and 0.20. Fig.[5]shows the 2D projection of momen- 




FIG. 8: SC order parameter, <I> versus hole concentration, 
S for different values of t± with J± = 0.1 The numbers in 
parentheses shown in figure are the values of {t±, J±). 

varies significantly and differently in the optimally doped 
and the overdoped regimes with the interplanar hopping 
t±. In the underdoped region, the effect of t± is small. 



FIG. 9: (color online) Plots of n(k) projected on a two- 
dimensional {kx, ky) plane where white (yellow) and gray 
(magenta) colors represent the highest and lowest spectral 
weights respectively. The left panel shows n{k) for fcz = 
and the right panel corresponds to k^ — tt, for t± = 0.05 (top) 
and 0.20 (bottom). The hole concentration here is 5 = 0.32. 



turn distribution function, 



n(k) with n(k) = (cj^^Ck^) 
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for 6 = 0.32. At this hole dopmg value there is a large 
variation in SC order parameter for the two values of 
t± (actually SC order is vanishing small a,t S = 0.32 for 
t± = 0.05) It is seen from the plots that there a sig- 
nificant variation of n(k) weights in the two momentum 
planes, viz. the = and kz = tt, as t± is varied. The 
spectral weight shifts from the kz = tt plane to the kz — 
plane as t± is increased from 0.05 to 0.20. Thus for larger 
t± , occupation of the pairs with fcz = is higher in the 
overdoped region. This favours the stability of SC phase 
as the kz = pairs are planar and contribute in develop- 
ment of the SC order. The kz = tt pairs are interplanar 
and hence are not expected to be key players as the SC 
correlations are essentially planar. We believe that this 
transfer of weight from kz — tt to kz — for larger tj_ 
helps in stabilizing the SC order in the overdoped region 
even though the disrupting effect due to larger interpla- 
nar hopping persists. In the optimally doped region, the 
momentum distribution profiles for different values of tj_ 
are found to be almost identical (plots not shown here) 
and hence the only effect of t± is to reduce the SC cor- 
relations. 

Regarding the finite size effects in our results, we would 
like to mention here that the magnitudes of various quan- 
tities calculated do show some dependence on the size of 
the lattice. The dependence of energy on lattice size is 
elucidated in details in Ref. [2^ However the main fea- 
tures of the key results here, e.g. the coexistence of AF 
and SC order, the effect of interlayer parameters on the 
properties etc. will remain qualitatively same with lattice 
size. 



IV. CONCLUSION 

We summarise our main results obtained using varia- 
tional calculations for a t-J bilayer as follows - coexis- 



tence of AF-SC is found to be more stable than the pure 
(d-wave) SC state at low values of doping (S < 0.14). 
Beyond this, of course, the SC state is found to have 
the lowest energy and remains stable upto a hole con- 
centration that is more than that obtained for the two- 
dimensional square lattice. Further, a detailed analysis 
of magnetic and superconducting properties yields the 
coexisting phase, not only energetically stable, but also 
supports a stronger AF and SC correlations. The third 
and possibly the most important result emerges when the 
effects of interplanar coupling parameters are invoked for 
discussion and it appears in the form of an asymmetry in 
the optimally and the overdoped region where the planar 
SC correlations are found to be more stable for a smaller 
and larger interplanar hopping respectively. However the 
interplanar exchange does not play a decisive role in SC 
correlations. 

Some of the other results obtained by us include a com- 
parison of the nature of magnetic correlations between 
the pure SC state and the AF-SC state. It is found that 
strong planar and interplanar AF correlations exists in 
the AF-SC state at and slightly away from half-filling, 
whereas in the pure SC state the magnetic correlations 
are very weak. We have also discussed the effects of inter- 
planar coupling on magnetic correlations in the AF-SC 
state where larger tj_ reduces planar magnetic correla- 
tions whereas larger J± enhances interplanar magnetic 
correlations. 
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